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Abstract.  A  simple  analytical  model  is  proposed  for  estimating  grain  boundary  mobility  during 
dynamic  recrystallization  in  metallic  alloys.  The  combined  effects  of  solutes  (solute  drag)  and 
second  phase  particles  (Zener  pinning)  on  mobility  are  considered.  The  approach  is  based  on  (and  is 
consistent  with)  a  recently  published  mesoscale  model  of  discontinuous  dynamic  recrystallization. 
The  dependence  of  grain  boundary  mobility  on  solute  concentration  and  particle  size  is  summarized 
in  the  fonn  of  two-dimensional  maps. 

Introduction 

Grain  boundary  migration  plays  an  important  role  in  dynamic  recrystallization  because  it  is  one  of 
the  main  parameters  controlling  the  final  grain  size  after  hot  working  of  a  material.  In  metallic 
alloys,  the  grain  boundary  mobility  M  is  generally  decreased  by  both  solute  atoms  (solute  drag 
effect)  and  second  phase  particles  (Zener  pinning)  [1],  Within  the  range  of  large  grain-boundary 
velocities  involved  in  classical  "discontinuous"  dynamic  recrystallization  (DDRX),  M  is  expected 
to  drop  monotonically  with  increasing  solute  concentration.  By  contrast,  the  local  grain  boundary 
migration  rate  either  falls  to  zero  or  remains  unchanged  according  to  whether  the  driving  force  is 
less  or  greater  than  the  Zener  pressure  Pz.  Such  discontinuous  behaviour  precludes  any 
straightforward  introduction  of  Zener  pinning  into  models  of  DDRX,  and  it  is  therefore  necessary  to 
estimate  an  average  value  of  the  mobility  pertaining  to  a  whole  grain  over  its  lifetime. 

Mesoscale  Power-Law  DDRX  Model 

In  a  mesoscale  "average  field"  model  for  DDRX  recently  published  [2,  3],  each  grain  is  considered 
as  a  spherical  “inclusion”  of  current  diameter  D  and  (homogeneous)  dislocation  density  p 
interacting  with  a  matrix  of  average  dislocation  density  p  .  The  grain  size  changes  according  to: 

dD/ds  =  (2Mx/s)(p- p),  (1) 


? 

in  which  s  and  s  are  the  strain  and  strain  rate,  respectively,  and  i  ~  pb  is  the  line  energy  of 
dislocations.  According  to  Eq.  1,  each  grain  grows  whenever  its  dislocation  density  is  less  than  p, 
and  shrinks  in  the  opposite  case.  The  evolution  of  dislocation  density  reflects  strain  hardening  and 
dynamic  recovery.  Various  equations  have  been  proposed  to  describe  such  phenomena;  most  of 
them  involve  numerical  computations.  In  the  present  work,  however,  the  following  simple  power 
law  is  used,  because  it  leads  to  closed-form  analytical  expressions: 


1 


dp/ds  =  Hv+1/ pv  , 


(2) 


in  which  v  >  0 ,  and  H  has  the  dimension  of  dislocation  density  (pm-2).  Note  that  Eq.  2  is 
associated  with  the  classical  macroscopic  (power  law)  Swift  equation  a  =  K(s  +  Sj  )n ,  in  which  K, 
s1?  and  the  macroscopic  strain-hardening  parameter  n  =  1/(2+  2v)  are  constants.  During  the  steady 
state  of  DDRX,  p  remains  constant.  Integration  and  combination  of  Eqs.  1  and  2  then  lead  to: 


D  = 


2Mt 


fpWV  k 

,H, 


yv  +  l  v  +  2  j 


When  the  grain  vanishes,  its  dislocation  density  reaches  its  maximum  value 


Pend 


v  +  2  _ 


V  +  1 


(3) 


(4) 


obtained  by  setting  D  =  0  in  Eq.  3.  (The  other  solution,  p  =  0,  corresponds  to  the  nucleation  of  a 
new  grain.) 

Grain  Boundary  Mobility  during  DDRX  of  an  Alloy  Containing  Precipitates 

When  stable,  incoherent  precipitates  are  present  in  the  material,  they  exert  a  drag  (Zener)  pressure 
Pz  that  opposes  grain  boundary  motion.  Setting  pz  =  Pz/t  , 

-If  |  p  -  p  |  <  pz ,  the  migration  rate  is  zero  because  the  boundary  is  pinned  by  precipitates,  and 
therefore  M  =  0; 

-If  |  p-  p  |  >  pz  ,  the  boundary  moves  with  mobility  M ,  (the  mobility  in  the  absence  of  Zener  drag), 

possibly  affected  by  solute  atoms  (see  below),  as  illustrated  in  Figure  1.  (However,  M,  is  likely  to 
depend  also  on  the  Zener  pressure,  a  question  which  is  left  for  later  investigation.) 


Figure  1.  Schematic  dependence  of  mobility  with  respect  to  the  dislocation  density  in  the  grain 
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The  average  mobility  pertaining  to  the  grain  over  its  lifetime  (or,  equivalently,  pertaining  to  the 
whole  set  of  grains  under  the  DDRX  steady  state)  is  then  given  by: 


P  max  |*  P  max 

M  =  S1/S,  inwhichS^I  M(p)D(p)“pvdp  and  S=l  D(p)2pvdp 

Jo  J  ( 


(5) 


o 


In  Eq.  5,  the  average  of  M(p)  is  weighted  by  the  current  surface  of  the  grain.  The  additional  factor 
pv  occurs  because  of  the  change  of  variable  from  time  t  to  dislocation  density  p.  Two  cases  are 
distinguished  for  evaluating  the  integral  for  Sp 

-If  P+Pz  <Pend>  Pz  ^P /(V+1)> 


Sj  =Mj 


P  P“Pz  p  Pmax 

J  D(p)2pvdp+M,J  D(p)2pvdp 


(6a) 


P+Pz 


-If  p  +  pz  >pend,  i.e.,  pz  >p/(v  +  l), 


Sl=Ml 


1 


P-Pz 


D(p)2pvdp 


(6b) 


Closed  fonned  expressions  are  then  obtained  for  the  above  integrals.  Figure  2  shows  that  M  drops 
rapidly  with  increasing  Zener  pressure  Pz  for  given  deformation  conditions,  i.e.,  a  given  p. 

However,  the  influence  of  precipitates  decreases  for  low  values  of  v,  i.e.,  when  strain  hardening  is 
strong  (per  Eq.  2). 


Pz/P 


Figure  2.  Dependence  of  the  grain  boundary  mobility  M  on  Zener  pressure  pz  (solid  lines)  and  the 
corresponding  linear  approximations  at  low  values  of  pz  (broken  lines).  Numerical  values  of  the 
slope,  k7,  are  given  for  v  =  0,  1,  and  4,  i.e.,  for  n  =  0.5,  0.25,  and  0.1,  respectively 
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It  appears  that  the  linear  approximations  of  M  given  by  the  initial  slopes  of  the  curves  in 
Figure  2  are  quite  acceptable  for  M/M,  >  0.3 .  This  is  relevant  because  DDRX  requires  sufficient 

grain  boundary  mobility  for  it  to  occur.  For  pz/p«  1,  first  order  polynomial  expansion  of  the 
integrals  in  Eq.  6a  leads  to  a  relation  for  M/Mp 

M/Mj  =l-kz  (pz/p)  in  which  kz  =  (3v  +  3)(3v  +  4)(3v  +  5)(v  +  l)3v+3/(v  +  2)3v+5  (7) 

It  is  interesting  to  note  that  a  linear  equation  similar  to  the  above  has  been  recently  proposed  by 
Hutchinson  et  al.  [11- 


Overall  Expression  for  the  Grain  Boundary  Mobility 

The  influence  of  solutes  on  the  grain  boundary  mobility  can  be  described  by  the  classical  equation: 


Mi  = 


M0 

l  +  PCm 


(8) 


pertaining  to  the  range  of  large  migration  rates  occurring  during  DDRX  [4],  Here,  M0  is  the 
mobility  in  the  pure  metal,  Cm  is  the  atomic  concentration  (atoms/nr)  of  solute  in  the  matrix,  and  [] 
a  constant.  It  is  more  convenient  to  write  Eq.  8  in  a  slightly  modified  form: 


M,  = 


M0 

1  +  a(Cm/Cs) 


(8a) 


in  which  the  maximum  {saturation)  solute  concentration  Cs  (which  is  a  function  of  temperature,  and 
may  be  modified  under  deformation  conditions)  has  been  introduced,  and  a  is  a  non-dimensional 
constant.  Contrary  to  the  case  of  Zener  pinning,  it  is  not  possible  here  to  use  a  linear  expansion  of 
Eq.  8a,  since  Cm/Cs  is  expected  to  vary  over  the  whole  range  between  0  and  1,  or  even  exceed 
unity  (supersaturation). 


The  grain  boundary  mobility  in  solute-containing  metals  undergoing  DDRX  is  therefore 
obtained  by  combining  Eqs.  7  and  8a,  viz., 


M  =  M0 


l~kz(pz/p) 

1  +  a(Cm/Cs) 


(9) 


Particle  Size  and  Solute  Concentration  Dependence  of  Grain  Boundary  Mobility 

Assuming  spherical  particles,  the  Zener  pressure  is  given  by  [4]: 

Pz=(7r/2)ynd2  (10) 

in  which  y  is  the  surface  energy  of  the  precipitates,  n  is  the  number  of  precipitates  per  unit  volume, 
and  d  denotes  their  diameter.  In  addition,  alloy  atom  conservation  leads  to  the  following  equation: 

Cm  +  ( V 6)  d3  n  Cp  =  C  (11) 
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3 

in  which  Cp  is  the  concentration  (atoms/m  )  of  solute  in  each  precipitate,  and  C  is  the  overall  alloy 
concentration.  Combining  Eqs.  10  and  11  yields: 


Pz=3y 


c-cm 

Cpd 


Eq.  9  can  then  be  written  in  the  fonn: 


!  _  gk  C/C„-Cm/Cg 


M  =  M 


o 


1  +  aCm/Cs 


with  8  = 


3y 


Tpcp/Cs 


(12) 


(13) 


Here,  8  is  a  constant  (with  units  of  length)  for  given  deformation  conditions.  The  dependence  of 
M/M0  on  d  and  Cm  can  then  be  plotted  in  a  two-dimensional  diagram.  However,  an  additional 

condition  must  be  fulfdled,  viz.,  L/d  >  (L/d)c  >  1 ,  in  which  L  denotes  the  average  center-to-center 

distance  between  precipitates,  and  the  critical  ratio  (L/d)c  means  that  the  Zener  pinning  theory  is 

likely  to  hold  only  for  sufficiently  widely-spaced  precipitates.  Using  n  =  l/L3  ,  Eq.  1 1  can  then  be 
recast  in  the  form: 


d  t6C-CmJ 

which  leads  to  the  condition: 

Cm  >  C  n  Cp/Cs 
cs  cs  6  (L/d)c3 


(14) 


(15) 


to  be  fulfilled  by  the  solute  concentration  Cm. 

An  example  calculation  of  the  grain-boundary  mobility  as  a  function  of  material  and  alloy 
parameters  is  shown  in  Figure  3.  For  this  example,  Cp  =2.23  x  1028  atoms/m3  was  estimated  for 

the  case  of  niobium  in  hexagonal  Ni3Nb  precipitates,  and  it  was  assumed  that  C/Cs  =  0.5 
(“undersaturated”  alloy  concentration)  and  (L/d)c  =  1.5  .  Other  input  parameters  are  specified  in  the 

caption.  In  the  cross-hatched  area,  the  grain  boundary  mobility  is  zero.  M  values  are  meaningful 
only  for  L/d>(L/d)c  =  1.5,  according  to  the  above  analysis.  Furthermore,  the  diagram  exhibits  a 

critical  particle  diameter  dc.  For  d<dc,  the  mobility  increases  with  Cm  (because  the  particle  spacing 
becomes  larger,  thus  weakening  the  Zener  pressure).  For  d  >  dc,  the  mobility  decreases  with 
increasing  Cm  because  solute  effects  predominate.  For  d  =  dc,  grain  boundary  mobility  is 
independent  of  the  solute  concentration. 


Concluding  Remarks 

A  model  has  been  developed  to  estimate  the  combined  effect  of  solutes  and  precipitates  on  the  grain 
boundary  mobility  under  conditions  of  dynamic  recrystallization.  Although  the  various  material 
parameters  are  not  well  known,  it  nevertheless  indicates  general  trends  and  orders  of  magnitude  that 
can  be  used  as  inputs  for  modeling  DDRX  in  metallic  alloys.  Conversely,  a  recently-published, 
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mesoscale  model  of  DDRX  enables  the  mobility  to  be  determined  from  the  steady-state  flow  stress 
and  average  grain  size  measurements  [3].  It  is  thus  also  possible  to  estimate  unknown  material 
parameters  from  such  data. 


C/Cs=  0.5 


Figure  3.  Diagram  showing  the  dependence  of  the  grain  boundary  mobility  (M/M0)  on  the  particle 
diameter  and  solute  concentration  for  an  overall  alloy  concentration  C/Cs  =  0.5.  Other  parameters 

assumed  for  the  model  calculations  comprised  y  =  0.5  J/m2;  kz  =  1.88  (associated  with  v  =  0  and  a 
macroscopic  strain  hardening  exponent  of  0.5);  t  =  3.306  x  10-9  J/m  and  p  =  510.2  /pm2  (both 
estimated  for  pure  nickel  at  900°C);  Cp  =2.23  x  1028  atoms/m3  (for  niobium  in  hexagonal 

Ni3Nb);  Cs  =  0.853  x  1028  atoms/m3  (for  niobium  in  nickel  at  900  °C);  and  a  =  2. 
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